On Rice's formula for stationary multivariate 
piecewise smooth processes 



K. Borovkov* and G. Last t 



September 21, 2010 



Abstract 



Let X = {Xt : t > 0} be a stationary piecewise continuous Revalued process that 
moves between jumps along the integral curves of a given continuous vector field, 
S C M. d be a smooth surface. The aim of this paper is to derive a multivariate 
version of Rice's formula, relating the intensity of the point process of (localized) 
continuous crossings of S by X to the distribution of Xq. Our result is illustrated 
by examples relating to queueing networks and stress release network models. 
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1 Introduction 

The classical Rice's formula going back to [TU] gives the intensity u(u) of crossings (orig- 
inally, upcrossings) of a given level u by a univariate continuous stationary Gaussian 

process X t in terms of the joint distribution of (X t ,X[) = (X ,X Q ), the process' value 
and its derivative at a fixed time (provided that the derivative exists in some suitable 
sense, e.g. in mean quadratic): 



where p(-, •) is the joint density (X , Xq) which is assumed to exist. Later on the result has 
been extended to more general classes of different iable (in some suitable sense) stationary 
processes, covering not only the first moments but also higher order factorial moments 
of the numbers of crossing, and even to more general settings for continuous random 
processes and fields. The formula proved to be quite useful in a number of applied areas, 
including signal processing, reliability, sea waves and others. For detailed accounts of the 
history of results of this kind and further bibliography, the interested reader is referred 
to HI], M and Chapter 3 in 0. 
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The case of processes with jumps and degenerate finite-dimensional distributions drew 
much less attention, although, from the applications' viewpoint, it is scarcely less inter- 
esting than the one of continuous processes. However, the heuristics behind the formula 
based on "Kac's counting formula" giving the number of crossings of a level u by a C 1 - 
f unction / on [0, 1] as 



(under a couple of further technical assumptions and denoting by 1A the indicator function 
of the set A), seems to be applicable in that case as well, provided that the process jumps at 
finite intensity and is smooth between the jump times. We note that fll -2ft is a consequence 
of Federer's coarea theorem (see e.g. (7.4.15) in pQ) 



applied to the function g(s) :— l{\f(s) — u\ < 5}. 

An analogue of (II. ip for the intensity v c { u ) of continuous level crossings by general 
univariate piecewise deterministic Markov processes that has the form 

u c (u) = \n(u)\p(u), 

where //(•) is the drift coefficient of the process and p(-) the density of X Q , was established 
in [6] (see also Theorem 12.51 below; one should mention here an earlier paper [I] where the 
case of Poisson shot-noise processes was considered). The result was used in [6] to obtain 
the asymptotic behaviour of the point processes of high level crossings (i.e. as u — > oo) in 
a number of interesting and important for applications special cases. 

The proof in [6] relied on the Markov structure of the process and in fact did not 
assume the existence of the density p — its existence was part of the assertion of the main 
theorem there. The natural question on whether Rice's formula for piecewise smooth 
processes can be extended to the multivariate and non-Markovian cases remained open. 
In the present paper we give a positive answer to it. 

The paper is organized as follows. In Section 2, we describe the main class of processes 
we will be working with and then present the main result together with its proof. Section 3 
presents examples to illustrate our main result. 

2 The main result 

First we will describe the main model of multivariate random processes X = {X t : t > 0} 
dealt with in this paper. The two key elements of the model are a point process N = 
{N(B) : B G £>(1R + )} of jumps in our process X (here and in what follows, B(-) denotes 
the class of Borel subsets of •) and a vector field /i : D — > W 1 defined on an open domain 
D C M d and specifying the dynamics of X between the jumps. Note that we allow trivial 
jumps (of zero size) as well. 

We assume that the following assumptions are satisfied. 
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(A.l) N is a stationary simple counting process on K + , which has a finite intensity 
X N := EN((0, 1]) and is such that N(R+) = oo a.s. 

The latter implies that the process A" is locally finite and hence its points can be 
enumerated in the increasing order. We denote them by < T\ < T2 < . . . and set 
To := for convenience (this is not a point of N a.s.). 

(A.2) fi G C\D). 

This assumption implies that there exist continuous functions t- : M d — > (— 00, 0) 

and t + : M. d — > (0, 00) such that, for any x G D, there exists a unique C 1 -function 
q(x, ■) : (t-(x),t+(x)) —> D satisfying the integral equation 

q(x,t) = x+[ fi{q{x,s))ds, t G (t_(x), t + (x)) (2.1) 
Jo 

(Picard-Lindelof theorem, see e.g. p. 8 in [8]). Moreover, for any fixed x G M. d there is 
a neighbourhood of (x,0) G R d+1 in which q(-, ■) will also be continuously different iable 
(Peano's theorem on dependence on initial conditions, see e.g. p. 95 in [S]). 

The integral curves q specify the dynamics of the process X between its jumps. 

(A. 3) Assume that, for any n > 0, one has Xx n G D, T n+ \ < T n + t + (XT n ) and 

X t = q(X Tn ,t - T n ), T n < t < T n+1 . (2.2) 

Moreover, X and are jointly stationary, i.e. the distribution of the bivariate process 
{(X s+t , N((s, s + t])) : t > 0} does not depend on s > 0. 

Next we will list assumptions involving the surface S of which the continuous crossings 
by X we are concerned with. The latter are defined as follows: we say that X has a 
continuous crossing of S at time s > if X s _ = X s G S and there is a 5 > such that 
X t ^S forte (s-5,s + 5)\{s}. 

(A. 4) Let S C D be the relative interior of a (d — l)-dimensional (not necessarily 
connected) C 1 -manifold with or without boundary, and {n(x) : x G S} be a continuous 
field of unit normals to S. Denoting by (•, •) the Euclidean scalar product in IR d , we 
assume that 

(n(x),/x(z)) 7^0, xeS. (2.3) 

Remark 2.1. Let r x := inf{t > : q(x,t) G S} be the first positive time the integral 
curve of \x leaving from x at time zero hits the surface S. It is not hard to see that, if S' 
is a compact subset of S, then from (A.2) and (A. 4) it follows that 

inf{r z : x G S'} > 0. (2.4) 

This (together with the fact that X jumps only finitely often in finite time intervals) 
implies that the times of continuous crossings of S through a compact subset of S cannot 
accumulate in finite time. 
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The times of continuous crossings of S by X form an at most countable set N c that 
will be identified with a random counting measure on [0, oo). Then 

$ C (C) := H(s,X s ) eC}= I l{(s,X s ) e C}N c (ds), C G B([0,oo) x R d ), (2.5) 

seN c 

defines a random (integer- valued) measure $ c on [0, oo) x M. d . For t > and 5" G B(M. d ), 
the random variable $ c ([0,t] x 5") need not be finite. However, if S' is a compact subset 
of S then (12 Ah implies that $ c ([0,t] x 5") < oo. Moreover, since N has a finite intensity, 

v c (B) := E$ c ((0, 1] x 5), fie £(M d ), (2.6) 

is finite, whenever B is a compact subset of S. Therefore z/ c (-) is a a- finite measure on 
B(WL d ). For any compact B C S, the point process $ c (- x B) is stationary. This is enough 
to derive the (refined) Campbell theorem stating that 

E J g(s, X s ) N c (ds) = J J g(s, x) ds u c {dx). (2.7) 

for any measurable function g : R + x M. d — > R + , cf. e.g. (1.2.19) in [3]. 

Remark 2.2. Assuming that A is a small enough open set to ensure that v c (SnA) < oo, 
observe that v c (S D A) _1 ^ c (-) can be interpreted as the distribution of the value of X at 
a typical time of continuous crossing of S" PI A. This is a particular instance of a Palm 
distribution, see e.g. [3j. 

The Palm measure ttq of the pairs of values of X just before and after a typical jump 
of X is defined by 
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(B) :=Ej2l{T n < l,X Tn _ ^ X Tn , (X rn _,X r J G 5}, 5 G £(M d x M d ). (2 A 



n=l 



Note that 7r (M d x R d ) < X N < oo. 

(A. 5) The distribution 7r of Xq has a continuous density p in a neighborhood of 5", 
and 

min{7r ((M d \ S) x S 1 ), tt (S x (R d \ S))} = 0. (2.9) 
Now we are ready to state our main result. 
Theorem 2.3. Under assumptions (A.1)-(A.5), one has the identity 

u c (B) = [ \(n(x),fi(x))\p(x)'H d - 1 (dx), B G B(R d ), (2.10) 
J SnB 

where H, d ~ l is the id — 1)- dimensional Hausdorff measure on lR d . 



Remark 2.4. As it will be seen from the first half of the proof of Theorem 12.31 the full 
continuity assumption on p (which is part of (A. 5)) can actually be somewhat weakened 
to the boundedness of p in a neighborhood of S and its right-continuity (in case the first 
number in (12. 9p vanishes) on S along the flow meaning that p(q(x, 0+)) = p(x), x G S. 



4 



In the one-dimensional case the above theorem simplifies to the following assertion. 

Theorem 2.5. In the case d = 1, assuming that S = {u} for some u G D such that 
ji{u) 7^ 0, and that (A.1)-(A.3) and (A. 5) are satisfied, one has 

Vc{[u}) = \fJi(u)\p{u). (2.11) 

Remark 2.6. In the Markovian case, representation ( 12. lip was established in [6]. More 
precisely, it was shown there that there exists a density p satisfying ( 12. lip . Due to the 
Markovian structure of the process, it was possible to derive the result under weaker 
technical assumptions. 

In the case when 

S = S u := {x G R d : x x = u} 

for some u G K, a continuous crossing of 5* is a continuous crossing of the level u by the 
first component of X. In this case Theorem 12.31 takes the following form. 

Theorem 2.7. Let assumptions (A.1)-(A.3) be satisfied andu G M. be such that S u C D 
and Hi(x) ^ for all x G S u , where [i\ is the first component of fi. Assume that (A. 5) 
holds with S = S u . Then, for B G B(R d ), 

u c (B) = J ■■■ J l B (u, x 2 , • • • , x d ) |/iiO, x 2 , . . . , x d ) | p(u, x 2 , . . ■ , x d ) dx 2 --- dx d . 

Remark 2.8. Theorem 12.71 is another, more straigthforward, generalization of ( 12. lip . 
Assume now that < v c (S u ) < oo and consider a "typical time" of a continuous crossing 
of the level u by the first component of X. Then the measure Q M (-) := i' c (S u )~ 1 b' c ({u} x •) 
describes the distribution of the other components of X at this time. This distribution 
can be interpreted in terms of the drift-modulated density p\ proportional to \[i\{x)\p[x) 
(assuming that E|/i!(X )| < oo). If (Yi, . . . , Y d ) is a random vector with density pi, then 
Q M is the conditional distribution of (Y 2 , . . . , Y d ) given that Y\ = u. 

Remark 2.9. Let k G {1, . . . , d} and assume that S = SxR d ~ k , where S C R k is a (& — 1)- 
dimensional smooth surface. Let {n(x) : x G S} be a continuous field of unit normals to 
S. Let X := (X«, . . . , X&) and Y := (X^ k+1 \ . . . , where X = (X^\ . . . , 

There is a one-to-one correspondence between the continuous crossings of S by the process 
X and continuous crossings of S by the process X. Equation fl 2 . 1 [) can be written as 

"c(B)= [ f\(fi(x),fi(x,y))\l B (x,y)p(x\y)n k ~ 1 (dx)F(Y edy), B G B(R h x Y), 

(2.12) 

where Y := M. d ~ k , jl(x, y) is the vector of the first k components of n(x, y), and x H- p{x\y) 
is the conditional density of X given that Y = y. In this form the result might be 
generalizable to other stationary pairs (X,Y). The process X should remain piecewise 
deterministic for given Y. But the process Y might take values in a more general space 
Y. In this paper we will make no attempt to establish such an extension of our results. 
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To prove Theorem 12.3} we will need an auxiliary result that requires some further 
notation. First of all, for our purposes it will suffice that that result would hold in a 
"local setting" , i.e. for S replaced with S l~l A, where A is a small enough open subset 
of IR d . As can easily be seen from the observation that we made after stating assumption 
(A. 2) and from (A. 4), if we understand by S such a "small piece" of the original surface, 
then the following will be satisfied: 

(A. 6) The surface S is connected and relatively compact, (12 .4p holds with S' = S 
and v c := u c (S) < oo. Furthermore, there exists a u > such that t + (x) > Uq for all 
x G S and, for any u G [0, u ], 

S u := {q(x,u) : x G S} 
is a C 1 -surface with a continuous field {n u (x) : x G S u } of unit normals to it satisfying 
mf{(n u (x),ii(x)) : x G S u , u G [0,u ]} > 0. (2.13) 
Moreover, it has a density p in a neighbourhood of S7o )UO ), where 

Sj := |J S u , I C R. 

Now denote by TV" the stationary point process of the times of all continuous crossings 
of S u by X. For any C G B([0, oo) x M d ), let $^(C) be the number of all s G iV c u such 
that (s,X s ) G C and := E$*([0, 1] x B), B G i3(M d ). 

Proposition 2.10. Under assumptions (A.1)-(A.4) and (A. 6), for any measurable 
function g : M. d — > R + , one has 

g{x)v u c {dx)= ( \(n(x),fi(x))\g(x)p(x)H d - l (dx) (2.14) 
j s u 

for ^-almost all u G [0, Uq\. 

PROOF. For any j > set Tj := Tj A 1 and, in particular, Tq := 0. For j > 1 we define 

Fix a £? G £>(IR d ) and assume that w G (0, u ). By definition, $"(/,- x B) > if and only 
if LjCiSuDB ^ 0. On the other hand, (12TT1 implies that x5)<l, so that 

x B) = i{i i n5 tt nB^}. 

Therefore 

oo 

$^((o, i)xB) = J2 HLj n s u n b ^ 0} 

and, for any w G (0, u ), 

/ $"((0,1) x B)du = V / n 5 U D B ^ ®}du. (2.15) 

Jo j =1 Jo 
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Now set 

Jj{v) := {t £ Ij : X t £ S [0tV) }, Uj{v) := {u £ (0,v) : L,- fl S u ^ 0}. 

Clearly, the last two sets are either simultaneously empty or are open intervals of the 
same length; in the latter case, put Uj(v) := inf Uj(v). Therefore, 



f 1{L 3 n S u f)B^ ®}du = [ l{q(X Uj{vh u - Uj (v)) £ B}du 

JO JUj(v) 

= [ l{X t EB}dt= [ l{X t eS {0iV) nB}dt, (2.16) 
Jjj(v) Jij(y) 

so that (12.15j) becomes 

/ $^((o, i)xB)du= / i{x t e S(p, v) n B}dt. 

Jo Jo 
Taking expectations on both sides of the last relation and using Fubini's theorem and 
stationarity of X, we obtain that 

fV pi 

/ v u c {B)du = E / i{x t e S(p,v) n B}dt = ¥(X e S {0 , v) n B). 

Jo Jo 

As functions of B £ B(M. d ), both sides specify a measure, and so the standard argument 
shows that, for any measurable function g : M. d — > M + , 



du J g(x)v™(dx) = / g(x)p(x)n d (dx). (2.17) 

Now we can assume without loss of generality that S admits a C 1 -parametrization 
(w\, . . . , it?d_i) z(wi, . . . , Wd-i), where (wi, . . . , Wd-i) varies in an open set W C 
For (wx, . . . , Wd-i) £ W and tt £ [0, uq], define 

^(wt, w d -i, u) := q(z(w!, w d -i),u), 

which, for a fixed u £ [0,Mo], will be a C 1 -parametrization of the "parallel" surface S u . 

Next we denote by <9j the operator of partial differentiation with respect to Wi, i = 
1, . . . ,d — 1, and let d^ip '■= d/du. A simple linear algebra calculations shows that the 
Jacobian Jip of if) = if>(wi, . . . , Wd-i, u) satisfies 

\JiP\ = \(n u (if>),d d if>)\H=\(n u {if>),»(if)))\H, 

where H 2 = H 2 (w\, . . . , Wd-i, u) is the determinant of the matrix [{diif>, djif))^ . . =1 d _ y 
However, for any fixed u £ [0, u ], H(wi, . . . , w d -i, u) dw\ ■ ■ ■ dw d -i is the surface element 
of S u in the coordinates (wi, . . . ,Wd-i), so that changing coordinates on the right-hand 
side of fl2TT71) yields 

/ g{x)p{x)l-i d {dx) = I g(if))p(if))\Jif)\dwi- ■ ■ dw d -\du 

g{if>)p{if>) \ (n u (if)), /i(if>))\ H dwi ■ --dw^i du 

Wx(0,v) 

du ( g(x)p(x)\(n u (x),fi(x))\n d - 1 (dx), (2.18) 

J S u 

which immediately implies the assertion of Proposition 12.101 □ 
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Remark 2.11. Assume that / is a real- valued C 1 -function defined on an open domain 
D C R d , with non- vanishing gradient and such that S u = {x G D : f(x) = u} for all small 
enough u. Such a function exists, at least for suitably small pieces of S. We may then 
apply Federer's coarea theorem (see e.g. (7.4.15) in PQ) on each open interval {T'^^T'^) 
to the level sets of the function t \-> f(X t ). While this would provide an alternative way 
for deriving ( 12.161) . we have preferred to give a direct argument presented in the above 
proposition. In a quite similar spirit the coarea theorem can be used to derive Rice's 
formula for smooth processes, see Section 11.4 in [T]. It was actually U. Zahle who first 
used in [12] the coarea theorem to prove Rice's formula for certain continuous processes. 
We also note in passing that the coarea formula could be used to establish (I2.18P as well. 
However, our more explicit argument yields additional information that is needed in the 
proof of Theorem 12.31 

Proof of Theorem 12.31 Since both sides of ( I2.10p are a-additive in B, it is no 
restriction of generality to assume that assumption (A. 6) is satisfied. Moreover, we can 
assume that S admits a smooth parametrization as in the proof of Proposition ^. 101 This 
is due to the fact that the surface S can be represented as a "mosaic" of "small pieces" 
for which the assumption will be satisfied owing to assumptions (A.1)-(A.4) on the 
original S. 

Furthermore, it is not hard to see that, to prove the theorem, it suffices to demonstrate 
that (I2.14p holds at u = for continuous and bounded g. We will show that by proving 
that, under the assumption 7ro((M d \S') x S) = 0, both sides of (12.141) are right-continuous 
at u = 0, as Proposition 12. 101 will imply then the desired result. The case when only the 
second term on the left-hand side of (12.91) turns into zero (i.e. no(S x (M. d \ S)) = 0) can 
be dealt with in exactly the same way by establishing the left-continuity of both sides 
of (12. 14p at u = in this situation (essentially via a time- reversal argument). 

Using the notation from the proof of Proposition 12.101 and setting 

h u (w) := (rau(^),AiWO)0WOpWOj w = ( W U ■ • ■ > w d-i), ip = ip(w,u), 
we have, for u G [0,wo], 

/ (n u {x)^{x))g{x)p{x)H d - l {dx) = j h u {w) H (w , u) U d ~ x (dw) . (2.19) 
j s u J 

As noted after stating assumption (A. 2), one has q G C 1 , and so n u (i/}(y,u)) is a contin- 
uous function of u, leading to 

h 0+ (w) = (n(z(w)),fi(z(w)))g(z(w))p(z(w)). 

Similarly, as u — > 0+, H 2 (w,u) converges to the value of the determinant of the matrix 
y(diz(w), djz(w))) i:j=1 d _ r Now the dominated convergence theorem implies that (12.191) 
converges to the right-hand side of (I2.14p at u — 0. 

To establish the desired right-continuity of the left-hand side of (I2.14p we assume that 
(I2.13P holds. Introduce the following point process $d on M + x M d x R d : 

oo 

$,(•) := X>{*t„- ^ X Tn }l{(T n ,X Tn -,X T J e ■}. 

n=l 
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Let u G [0, Mo] and t > 0. A continuous crossing of S u can only occur on a trajectory of X 
that arrives at the surface from the inside of 5(o, u ) along an integral curve of [i (cf. (I2.13P ). 
Therefore each such crossing of S u should be preceded by an entry to Sp.u), either along a 
drift line or by jump. Taking into account the possibility of having X e S[o lU ), we obtain 
the bound 

N:([u,t + u]) < N c ([0,t + u}) + $ d ([0,t + u] x (R d \S [oM ) x S [0)U] ) + 1. 
Therefore 

+ g(X s )N?(ds) < f + (g(X s )+e(u))N c (ds) 
Jo 

+ g*[$ d ([0,t + u] x (M. d \S [0M ) x S [0>u] ) + 1], 

where g* : = sup x . g(x) and 

e(it) := sup{\g(q(x,v)) - : x E S, < v < u] -> as m-^0+ (2.20) 

due to the uniform continuity of the mapping (x, u) H- g(q(x, u)) on 5 x [0, u ], S denoting 
the closure of S. 

Now taking expectations on both sides of the obtained inequality and using Campbell's 
formula (12. 71) yields 



t J g{x)v u c {dx)<{t + u) J (g(x)+e(u))v c (dx)+g*((t + u)n ((R d \S [0M ) x S [0 , u] ) + l), 

where we also used Campbell's theorem for vr (-) = E$ d ([0, 1] x ■). After dividing by t 
and letting t — > oo, we obtain 



g(x)v u c (dx) < J g(x)v c (dx) + e(u)u c (S) + g*n ((^ \ S M ) x S [0>u] ). 

In view of (12.201) and the fact that the assumption 7r ((lR d \ S) x S) = implies that 
7r ((M d \ S[o )U ]) x S[o >u ]) — > as u — > 0+, this leads to 



limsup / g{x)v™(dx) < / g(x)v c (dx). 

To derive the converse inequality, we start with the observation that any continuous 
crossing of S in [0, t] is followed either by a continuous crossing of S u or by a jump from 
S(o,u] t° its complement within the time interval [0, t + u], so that 

NZ([0,t + u}) >N c ([0,t])-$ d ([0,t + u] xS (Q , u] x (R d \S {0 , u] )) - 1. 

Next, similarly to our argument above, we obtain 

s(xK(cfe) > y g{x)v c {dx) - e(u)u c (S) - g*ir Q (S { o !U ] x (M d \ 5 (0 , u] )). 

Since lim u ^o+ S(p jU ] = 0, it is clear that the continuity of 7Tq implies now that 



liminf / g(x)u^(dx) > / g(x)v c (dx), 

u->0+ J J 

which completes the proof of the theorem. □ 
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3 Examples 



In this section we will present two examples showing possible uses of our main result. 

Example 3.1. Consider a general queueing network model with d servers operating in 
stationary regime, with arrivals of customers (possibly in batches) to the network being 
governed by a stationary simple point process. Each customer, upon completion of its 
service at node j G {1, . . . , d} of the network, proceeds to another node for further service 
or leaves the network, according to some routing mechanism. All the arrival, transition 
and departure times form a stationary point process N, and it is at these times that the 
state of the process X t = (X t , . . . , X[ ) G lR d describing the residual workloads on the 
nodes can change by a jump. Between the events, the values of X t decrease according to 
the relation -^X t = Li(X t ) for some C 1 -function li : !BL d — > Wt, so that the service rate at 
node j can depend on the residual workload at the node and, moreover, it can even depend 
on the workloads at other nodes i ^ j as well. To make this description compatible with 
the assumptions in Section 2, we allow X^ < interpreting as the residual workload at 
node j at time t the value max{X^\ 0}, and let D := M. d . 

For i G {1, . . . , d} let Si := {x = (xi, . . . , Xd) G M. d : X{ = 0}. Then the continuous 
crossing of the surface Si corresponds to server i becoming idle. Let Vi[B) denote the 
intensity of these crossings through a point in B G Provided that the assumptions 

of Theorem 12.71 are satisfied, we obtain 



for B G B(M d ), where x % := (x±, . . . , 0, . . . , Xd) and fa is the ith component of 
li. The normalization of (13.11) yields the (Palm) distribution of the network at a typcial 
departure time from node %. 

Note that assumptions (A.1)-(A.3) are rather mild and that Theorem 12.71 also re- 
quires Hi{x) < for x G Si. In assumption (12. 9p only the condition no(Si x (M d \ Si)) = 
is of relevance. This assumption says that if there is a jump at an instant when server i 
becomes empty, then the workload of this server is not allowed to increase by this jump, 
neither by an internal transition (including feedback) nor by an external arrival. Again, 
this is a rather weak assumption. 

We can also consider the "composite surface" S := [JiS^, where S^) is the set of all 
x G M d with Xi = and Xj ^ for j ^ i. (Under weak assumptions any continuous cross- 
ing of Si is also a continuous crossing of >%)•) Theorem 12 . 71 provides the Palm distribution 
of the residual workloads at the time when one of the servers becomes idle while all the 
others are still working. The probability of server i becoming idle given a typical instant 
when (exactly) one of the servers becomes idle is then given by v c (S(i))/is c (S). We skip 
further details. 

Example 3.2. The classical stress release model in seismology (see e.g. [7] and references 
to earlier work therein) is a piecewise deterministic Markov process X t representing the 
level of "stress" at a seismic fault at time t. The value X t continuously increases at a 
linear rate due to the tectonic loading of the fault and drops by random jumps when 
the stress discharges by way of earthquakes that occur at random times whose intensity 
is given by ip(X t ) for some suitably chosen increasing risk function (e.g. ip(x) = e^ x 



MB) 



1 



1 



J-b(^) \fJ>i(x l ) | p{x l ) dxf ■ dxi-xdxi+i ■ ■ ■ dx d , 



(3.1) 



10 



for some /3 > 0). Note that the remote measuring of stress levels at seismic faults is an 
extremely difficult problem, so the value X t is usually not observable. All the information 
on the process one can have access to is contained in the times, locations and magnitudes 
of jumps. 

A more interesting multinode analog of the model was discussed in [5], where it was 
demonstrated, in particular, that already a two-node stress release network can reproduce 
the famous Omori's law for the intensity of earthquake aftershocks. 

In the multinode model, the values of the components of the random process X t = 
(X^\ . . . ,X^) G M. d represent the time t stress levels at individual seismic faults j G 
{1, . . . , d} constituting a local fault system. Between jumps, the dynamics of the process 
are given by 4zX t = ji for a constant vector u G M. d , Note that one can have jij < which 
corresponds to tectonic unloading of stress at node j (of course, we can consider a more 
general model with a variable u as well; similar remarks apply to all the other elements 
of the model construction). Jumps ("seismic events") occurrence at node j is driven by 
a Markovian random mechanism with the probability of a jump occurring at the node in 
the infinitesimal time interval dt given by ipj(X t _)dt for a given risk function ipj(x). 

When the nth seismic event occurs at node j (say, at time l} jn ), the value of stress at 
the node changes by a random quantity n = 1,2,..., which may be assumed to be 
i.i.d. random variables. Moreover, the stress levels at other nodes can also change at that 
instance: for a given constant transfer matrix (r^) G R dxd , one has Xj} = X£? _+r ii £ 7> , 
i 7^ j (for more detail, see [5]). 

One of the main problems one hopes to be able to solve in mathematical seismology 
is to give advanced earthquake warnings. Within the framework of the multinode stress 
release model, that warning would have to be given at the time when the cumulative 
jump intensity £)*=i V'P^) exceeds a given threshold u > 0. That is, we are looking at 

continuous crossings of the surface S := {x : Ylj=i ^( x j) = u } by our process X t . Our 
main result allows to find the distribution of X t at the (typical) time of such crossing and 
hence, for example, to derive the probability for a given fault to trigger the forthcoming 
seismic event. 
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